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Abstract 

We consider iid Brownian motions, Bj{t), where Bj{0) has a rapidly decreasing, 
smooth density function /. The empirical quantiles, or pointwise order statistics, 
are denoted by Bj-n{t), and we consider a sequence Qn{t) = Bj^j^y^{t)^ where 
j{n)/n — )• a G (0, 1). This sequence converges in probability to the a-quantile of 
the law of Bj{t). We first show convergence in law in C[0,oo) of Fn = n^^'^iQn — q)- 
We then investigate properties of the limit process F, including its local covariance 
structure, and Holder-continuity and variations of its sample paths. In particular, we 
find that F has the same local properties as fBm with Hurst parameter H = 1 /4. 

AMS subject classifications: Primary 60F05; secondary 60F17, 60G15, 60G17, 
60G18, 60J65 

Keywords and phrases: quantile process, order statistics, fluctuations, weak 
convergence, fractional Brownian motion, quartic variation 

1 Introduction 

The classical example of a Brownian motion occurring in nature is the pollen grain on the 
surface of the water. Bombarded by the much smaller water molecules that surround it, the 
pollen grain embarks on a random walk with numerous, tiny steps. Applying a law of large 
numbers type scaling to its path leaves it sitting right where it is, since its trajectory has 
mean zero. But applying a central limit theorem type scaling, which exposes its fluctuations, 
reveals the Brownian path that we observe in nature. 

If we observe a large collection of pollen grains, and approximate their density with a 
continuous function, then we might expect this density function to evolve according to the 
diffusion equation. In this setting, the diffusion equation can be derived by assuming that 
each individual pollen grain moves as an independent Brownian motion. 
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However, if the pollen grains are close to one another, as they should be if we are 
approximating their density with a continuous function, then their motions are certainly 
not independent. They are interacting with one another through collisions. Moreover, the 
collisions would evidently provide each individual pollen grain with a drift, pushing it toward 
regions where the pollen density is lower. In other words, the individual pollen grains would 
not be performing Brownian motions. The question then arises, what do the trajectories of 
these colliding pollen grains look like on the mesoscopic scale (that is, under central limit 
theorem type scaling)? 

In this paper, we consider the simpler, one-dimensional question of the scaling limit of 
colliding Brownian motions on the real line. To motivate our model of the collision process, 
first consider two physical particles with equal mass and velocities Vi and V2- If these 
particle have an elastic collision, then they will exchange velocities. The effect of this will 
be that they exchange trajectories at the moment of collision. If Xi{t) and X2{t) describe 
their trajectories without any interaction, then their trajectories in the presence of collisions 
will be min{xi{t) , X2{t)} and ma.x{xi{t),X2{t)}. If we extend this reasoning to n particles, 
then we are led naturally to consider the empirical quantiles, or order statistics, of the non- 
interacting trajectories, xi{t), . . . ,a:„(t). Two concepts, therefore, which are central to our 
approach, are the notions of an a-quantile of a probability measure, and the order statistics 
of a family of random variables. 

Let z/ be a probability measure on M with distribution function $^(x) = i/((— oo,x]). 
Given a G (0, 1), an a-quantile of z/ is any number q such that $i,(g— ) < a < ^u{(l)- It is 
easy to see that u has at least one a-quantile, given by g = inf{x : a < ^^{x)}. It is also 
clear that if $1, is continuous, then $,y(g) = a for any a-quantile, q. 

Given random variables, let a be a (random) permutation of the set 

{1, . . . , n} such that -^cr(i) < ■ ■ ■ < ^cr(n) a.s. For 1 < j < n, we define the j-th order statistic 
of X = (Xi, . . . ,Xn) to be X„Q), and denote it by Xj,n. Note that —Xj,n = {—X)(^n-j+i):n- 

Now fix a G (0, 1) and let i? be a one- dimensional Brownian motion with a random initial 
position. For simplicity in the calculations that are to come, we shall assume that -B(O) has 
a density function / that is a Schwartz function. That is, / G C°°(M) and 

sup(i + |xr)i/('")(x)i < 00, (1.1) 

for all nonnegative integers n and m. We will also assume that /(x) dx has a unique a- 
quantile, ^(O), such that f{q{0)) > 0. 

Let {Bj} be an iid sequence of copies of B. For fixed n, the trajectories Bi, . . . ,Bn denote 
the paths of a system of n particles with no interactions. When these particles are allowed 
to interact through collisions, their trajectories are given by the pointwise order statistics, or 
empirical quantile processes, Bi-n, ■ ■ ■ , Bnm- More specifically, Bj^n is the process such that, 
for each t > 0, Bj.n{t) is the j-th order statistic of {Bi[t), . . . ,Bn{t)). Note that Bj-n is a 
continuous process. We shall fix a sequence of integers {j{n)}'^^i such that 1 < j(n) < n and 
j{n)/n = a + o(n~^/^). We then consider the sequence of empirical quantile processes {Qn} 
given by Qn = Bj(^n) n- Our first result concerns the convergence in C[0, 00) of the sequence 
{Qn}- To this end, we begin by defining the (deterministic) limit process, and show that it 
is continuous. 

Let u{x,t) denote the density of B{t). It is well-known that our assumptions on / are 
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sufficient to ensure that u E C°°(M x (0, oo)), diu E C(R x [0, oo)), and diu{x, 0) = /(^^(x) 
for all j > 0. Moreover, diu{x,t) = E^[f^^\B(t))] and u satisfies the diffusion equation, 
dtu = {l/2)dlu. For each t > 0, the function u{-,t) is strictly positive, which implies 
that u{x,t)dx has a unique a-quantile, q{t). The following lemma is easily derived by 
differentiating the defining equation for q{t)] its proof is given in the appendix. 

Lemma 1.1. The function q is in C[0, oo) fl C°°(0, oo) and satisfies 

d^u{q{t),t) 

for all t > 0. 

Remark 1.2. Let C''[0, oo) denote the space of functions g : [0, oo) — ?■ M such that g^^'^ has a 
continuous extension to [0, oo), for all j < k. Also, let C°°[0, oo) = nk>oC''[0, oo). It is easy 
to see from Lemma [LI I that since u{q{t),t) > for all t > 0, and dlu E C(R x [0,oo)), it 
follows that q E C°°[0, oo). That is, q^''^ has a continuous extension to [0, oo) for all k. 

Remark 1.3. The diffusive flux at time t, in the positive spatial direction, across an 
arbitrary moving boundary s(t) is given by — (1/2) dxu{s(t),t) — u{s(t),t)s'(t). The ffist 
term corresponds to Pick's ffist law of diffusion, and the second term comes from the motion 
of the boundary. By (11. 2p . it follows that q{t) is the unique trajectory starting at g(0), across 
which there is no diffusive flux. The vanishing of the flux can also be seen by noting that 

lit) 

u{x, t) dx = a, 

-oo 

which follows from the definition of q{t). 

It will follow from later results that Qn — ?■ g in probability in C[0,oo). Our primary 
interest, however, is with the fluctuations, -F„ = n^^'^{Qn — q)- Our objective in this paper is 
twofold: (i) to establish the convergence in law in the space C[0, oo) of the sequence Fn to a 
limit process F, and (ii) to investigate properties of the limit process F, including its local 
covariance structure, and the Holder continuity and variations of its sample paths. 

Regarding (i), we establish the following result, whose proof can be found in Section Ei 
(The notation X„ ^ X means X„ — > X in law.) 

Theorem 1.4. There exists a continuous, centered Gaussian process F with covariance 
function 

^^^^^^_ P{B{s)<q{s),B{t)<q{t))-a^ ^ (1.3) 
u{q{s),s)u{q{t),t) 

such that Fn ^ F in C[0, oo). 

It is worth pointing out here that the limiting process F is not deterministic at time 
t = 0. In fact, F|F(0)p = {a - a^)f{q{0)y\ 

Regarding (ii), we derive, in Section [31 several key properties of the limit process F. We 
will show that, on compact time intervals, E\F{t) — F{s)\'^ is bounded above and below by 
constant multiples of |t — s|^/^. In particular, the paths of F are almost surely locally Holder 
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continuous with exponent 7 for all 7 < 1/4. (See Corollaries 13.21 and Corollary 13.31 ) We will 
show that F is locally anti-persistent. More specifically, nearby, small increments of size At 
have a negative covariance whose order of magnitude is —\t — s\~^^'^At^. (See Corollaries 
13.41 and 13.51 ) We will also show that F has a nontrivial, deterministic quartic variation, 
given by (6/7r) |M(g(s), s)|~^ (is. (See Theorem 13.71 ) All of these are local properties 
that F shares with 5^/^, the fractional Brownian motion (fBm) with Hurst parameter 
H = 1/4. (Recall that fBm, 5^, is a centered Gaussian process with 5^(0) = and 
E\B^(t) — B^{s)\'^ = |t — sp^, where H G (0, 1).) These local properties are also shared 
with the solution to the one-dimensional stochastic heat equation, which was recently studied 
in [T7] and [1]. On the other hand, the global properties of F can be quite different from 
B^/^, which we will illustrate at the end of Section [3l via the special case where f{x) dx is a 
standard Gaussian distribution. 

A model similar to this was studied in [7j. That model consists of a countably infinite 
collection of real- valued stochastic processes Xi(t), i E Z, t > 0. The points {xi(0)}jgz form 
a unit Poisson process on the real line, conditioned to have a point at 0, and labeled so that 

• • ■ < x_2(0) < a:_i(0) < a:o(0) = < a:i(0) < ^2(0) < ■ ■ ■ 

The processes {xi{-) — Xj(0)} are independent, standard Brownian motions. This family of 
processes represents the motion of the particles without collisions. By counting upcrossings 
and downcrossings, the motion of a "tagged particle" in the collision system can be defined. 
Informally, the tagged particle y{t) is defined as follows. It begins with y{0) = a;o(0), and 
then continues as y{t) = Xo{t) until the first time the path of Xq intersects one of its neighbors, 
Xi or x„i. At this point, y adopts the trajectory of the neighboring particle, and follows this 
trajectory until the first time it meets one if its two neighbors, and so on. 

Of course, when two Brownian particles meet, they intersect infinitely often immediately, 
which makes it difficult to carry out the above informal construction. This is why upcrossings 
and downcrossings are used instead. (See also [5].) Note that, with only finitely many 
particles, the empirical quantiles serve as tagged particles, without any need for counting 
crossings. 

By the scaling property of Brownian motion, the process n~^y{nH) has the same law 
as a tagged particle in a system initially distributed according to a Poisson process with 
a density of n particles per unit length. This is analogous to our centered empirical 
quantile process Qn{t) — Qi^)- Multiplying by ra^/^ shows that our process F„ is analogous to 
yA{t) = y{At)/A^'^, where A = . In [7], it is shown that ?/a(1) = y{A)IA^I^ 
as A —7- 00, where is a standard normal random variable. This implies that for fixed t > 0, 

VAit) = t'/^'yAtil) t'/'{2/7rY/^N ^ {2/7rY/'B'/\t), 

where = denotes equality in law. In [5], a substantially stronger result was proven, one 
case of which shows that yA =^ (2/7r)^/'^i?^/^ in C[0, 00). The proof of tightness in [5] relies 
heavily on the special properties of the initial Poisson distribution, under which the resulting 
particle system has a stationary distribution. 

Similar results hold for the simple symmetric exclusion process (SSEP) on the integer 
lattice. It has been known since the work of Arratia [1] and Rost and Vares [H] in the 1980s 
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that the fluctuations of a tagged particle in a one-dimensional SSEP started in equilibrium 
converge, in the sense of finite-dimensional distributions (fdd), to fBm 5^/^. The proof of 
tightness, however, was established only recently by Peligrad and Sethuraman [13] in 2008. 
As with the proof in [5], the proof of tightness for the tagged particle in [13] also relies on 
the initial distribution of the particles and does not extend to the non-equilibrium case. 

The motion of a tagged particle in a non-equilibrium SSEP was studied in Jara and 
Landim [8]. Theorem 2.6 in [8] establishes a result on the fluctuations of a tagged particle in 
SSEP which is analogous to our Theorem II. 4[ except that in [^ convergence is established 
only in the sense of fdds. 

It should be noted that the proof methods in [5] do generalize to the non-equilibrium 
case when one considers the current process, rather than the tagged particle. In 2005, 
Seppalainen [15] studied the current process in a system of independent random walks in 
a non-equilibrium case. Tightness for this current process (for a certain restricted class of 
initial profiles) was proved by Kumar [TT] in 2008 by extending the proof of Proposition 5.7 
in [g. 

In [18] , we considered what is very nearly a special case of the present model. Suppose 
-B(O) = a.s. and j{n) = [{n -|- 1)/2J , so that Qn is the median process, a = 1/2, and q = 0. 
In that case, it was shown in [18] that F„ ^ X in C[0, oo), where X is a centered Gaussian 
process with covariance function 



The main difficulty in establishing this result was the proof of tightness. Unlike [5], we 
did not have the benefit of the Poisson initial distribution. Instead, we proved tightness 
by making direct estimates on P{\Qn{t) — Qn{s)\ > The method of estimation was 

essentially a four-step procedure. First, we established a connection between this probability 
and a certain random walk. Second, we derived estimates for probabilities associated with 
this random walk - estimates which depend on the parameters of the walk. Third, we 
estimated those parameters in terms of the dynamics of the original model. And fourth, we 
put all of this together to get separate estimates on P{\Qn{t) — Qn{s)\ > n'^^'^e), depending 
on whether n~^/'^e <^ |t — s|^/^, n~^/^e f« |t — s|-^/^, or n~^/'^e ^ \t — 

Regarding the convergence portion of the present paper, our proof of tightness is 
motivated by this rough outline. It is our hope that the new ideas and techniques developed 
here in relation to tightness, as well as the methods employed to study the properties of the 
limit process, will be applicable to more general colliding particle models. 

2 Finite-dimensional distributions 

We summarize here the definitions and assumptions from Section [1] Let 5 be a Brownian 
motion such that -B(O) has a density function / G C°° satisfying ( II. ip for all nonnegative 
integers n and m. Let {j(n)}5^]^ be a sequence of integers satisfying 1 < j{n) < n and 
j{n)/n = a + o(n~^/^), where a G (0,1). Assume the probability measure f{x)dx has a 
unique a-quantile, q{0), satisfying f{q{0)) > 0. 




(1.4) 
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Let {Bj} be a sequence of iid copies of B. We define F„ = n^^'^{Qn—q), where Qn = Bji^ny.n 
and q{t) is the unique a-quantile of the law of B{t). The fact that F„ G C[0,oo) follows 
from Lemma ILII 

In this section, we begin the proof of Theorem 11.41 by noting that the convergence of the 
fdds follows as an immediate corollary of a multi-dimensional quantile central limit theorem, 
which is stated below as Theorem 12.11 The proof of the quantile CLT is a straightforward 
exercise in the application of the Lindeberg-Feller theorem, and is given in the appendix. 

To state the quantile CLT, we first need some preliminaries. Let X be an M'^-valued 
random variable. It will be convenient to denote vector components with function notation, 
X = {X(l),...,X{d)). Let = P(X(j) < x) and = P{X{i) < x,X{j) < y). 

Fix a G (0,1)'^ and assume there exists g G M'^ such that = a{j) for all j. Also 

assume that, for all i and j, exists and is strictly positive, and that Gij is continuous 

at {q{i),q{j))- 

Let {Xn} be a sequence of iid copies of X. If k G {1, . . . then X^^-n G shall 

denote the componentwise order statistics of That is, is the K(j)-th 

order statistic of (Xi(j), . . . , 

Theorem 2.1. With the above notation and assumptions, define the matrix a G M'^^'^ by 

^ Gij{q{i),q{j)) -a{i)a{j) . . 

If K{n) = (/t(n, 1), . . . , /t(n, d)) G {l,...,n}'^ satisfies K{n)/n = a + o(n~^/^), then 
Tn}^'^{X^(^n):n — q) ^ N, where N is multinormal with mean zero and covariance matrix 
a. 

Corollary 2.2. There exists a centered Gaussian process F with covariance function p, given 
by (11.31) . such that Fn ^ F in the sense of finite- dimensional distributions. 

Proof. Given < ti < ■ ■ ■ < t^, apply Theorem O with X = {B{ti), . . . , B{td)) and 
K{n) = {j{n)J{n),...J{n)). □ 

Remark 2.3. Note that one of the assumptions in the above theorem is that exists 
and is strictly positive. In particular, must be continuous at q{j). That is, the theorem 
does not apply when the distribution of X{j) has a jump at the quantile point q{j)- This, 
however, is never the case for the random variables considered in this paper, since for any 
initial distribution, the solution of the heat equation is continuous for any positive time. 

Remark 2.4. To give some limited intuitive sense to (12. ip . let us consider what it implies 
for stochastic processes. Let : t G T}'^^-^^ be a sequence of iid, measurable, real- 

valued copies of a stochastic process X(t), parameterized by some general index set T. 
Let ^{x,t) = P{X{t) < x). Suppose that is continuous and strictly increasing, 

and, for a G (0,1), define q{a,t) so that $(g(a,t),t) = a. Let k : (0, 1) x N N 
satisfy n{a,n) G {l,...,n} and K{a,n) = an + o(n^/^) for each fixed a. Let X^-nit) 
denote the K-th order statistic of Xi(t), . . . ,X„(t), and define Q„(a,t) = Xf^(^a,n):n(t). Let 
F„(a, t) = n^/'^{Qnia, t) - q{a, t)). 

From Theorem 12. H we find that if we make the additional assumptions that dx^{x, t) > 
for all X such that < $(x,t) < 1, and that the functions {x,y) i— )■ P{X{s) < x,X(t) < y) 
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are continuous for each fixed s and t, then the fdds of the two-parameter processes Fn 
converge to those of a centered Gaussian process F on (0, 1) x T with covariance function 

ma, s)FW, t)] = (2-2) 
' ^ ' ^ " dM<l(a,s),s)dM<l(P,t),t) ^ ' 

In this paper, we are considering not only the special case in which the processes are Brownian 

motions, but also the special case in which a is fixed. If we fix t, however, then a familiar 

process emerges. Namely, for fixed t, the process a h- )■ F{a,t) has covariance function 

E[F{a, t)F{P, t)] = ^ , ^^rf ^ ^ • (2-3) 

In other words, dx^{q{-,t),t)F{-,t) is a Brownian bridge. This phenomenon stems from 
the relationship between the quantiles and the empirical processes. That is, let = 
iE%ikx,{t)<x}. Then 

n^/^($"(Q„(a,t),t) - $(g(a,t),t)) = n^/^{n-^K{a,n) - a) = o(l). 

Hence, under suitable assumptions on and its derivatives, we may use a Taylor 

polynomial for about q{a,t), obtaining 

ni/2($"(g„(«, t),t)- $(g„(«, t),t)) = -ni/2($(g„(a, t),t) - t), t)) + o(l) 

= -9,$ t), t) + t) ^0(1) + o(l). (2.4) 

For fixed t, it is well-known that r2^/^($"(-, t) - converges to i?°(<l>(-,t)), where B° 

is a Brownian bridge (see Billingsley [2], for example). Hence, the left-hand side of (12.41) 
converges to B°{^{q{a,t),t)) = B°{a), which explains the covariance function (12.31) . 

The connection to the empirical processes given by (12. 4p is still valid, of course, even when 
t is allowed to vary. In this case, (12. 4p and (12.21) show that the two-parameter fluctuation 
processes ?7,^/^($" — $)—)■ \^ in the fdd sense, where V is centered Gaussian with covariance 

E[V{x, s)V{y, t)] = P{X{s) < X, X{t) < y) - P{X{s) < x)P{X{t) < y). (2.5) 

At this point, it is interesting to compare this with the results in Martin-Lof [12]. When X 
is a Markov process, [12] considers the convergence, in the fdd sense, of fluctuations such 
as these, although for a slightly different model in which we have ~ Poisson(?T,) particles 
instead of n particles. This subtle difference produces slightly different fdds in the limit. 
(For example, in the fixed t case, the fluctuations of the empirical distributions of a Poisson 
number of variables converges to a Brownian motion rather than a Brownian bridge, which 
can be seen either from the heuristic demonstration below, or as a special case of Corollary 
2 in [12].) 

More specifically, let $"(x,t) = ^ Yl!j=i '^{Xj{t)<x}- Note that 

^i/2($n _ $) = {N/nf'^N^'\^^ - $) + n^'\N/n - 1)$. 

Hence, n^/^($" — $) — )■ in the fdd sense, where V = V + Z$, and Z is a standard normal 
random variable, independent of V. Using (12.51) . it follows that V is centered Gaussian with 
covariance 

E[Vix, s)Viy, t)] = PiXis) < X, X(t) < y). 
This, of course, is precisely the result we obtain from Corollary 2 in 
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3 Properties of the limit process 

Before proceeding to the proof of tightness, we first estabhsh some of the properties of the 
hmit process. In what follows, C shall denote a positive, finite constant, that may change 
value from line to line. 

Theorem 3.1. Recall p given by ( \1.3\\ . For each T > 0, there exist positive constants 
S, Ci, C2 such that whenever < s < t < T and \t — s\ < 6, we have 

(l) < dsp{s,t) < C2\t-s\~'/\ 

(11) -C2\t- s\~^l^ <dtp{s,t) <-Ci\t- s\-^l'^, and 
(m) -C2\t - s|-3/2 < dlp{s,t) < -Ci\t - 

Proof. Define F{t) = u{q(t),t)F(t) and note that F is a centered Gaussian process with 
covariance function 

Pis, t) = P{Bis) < qis), Bit) < qit)) - a\ (3.1) 

We will first prove the theorem for p instead of p. For later purposes, it will be convenient 
at this time to introduce the auxiliary function 

/q{s)+w /•q{t)+z 
/ p{t - s,x,y)dlu{x,s)dydx, (3.2) 
-00 J —00 

where i >0 and p{t,x,y) = (27rt)~-'^/^e~*^^'~^^^/^*. In this proof, we use only the case i = 0. 
Later, in Section [6l we will make use of the general case for arbitrary i. We now compute 
that 



/q{s)+w 
pit — s, X, q{t) + z)dl.u{x, s) dx 
-00 



-00 

/q{s)+w i-q{t)+z 
I dtp{t — s, x,y)dl.u{x, s) dy dx 
00 J —00 

/q{s)+w 
p{t — s, X, q{t) + z)dl.u{x, s) dx 
'OO 

i>q{t)+z /'q(s)+w 

52 „ 



(3.3) 



d^p{t — s, X, y)dl,u{x, s) dx dy. 



Using integration by parts, 

fq{s)+w 

dlp{t — s, X, y)dl.u{x, s) dx 
dxp{t - s, q{s) + w, y)dl.u{q{s) +w,s)- p{t - s, q{s) + w, y)d'^^u{q{s) + w, s) 

fq{s)+w 

pit — s, X, y)dl^'^u{x, s) dx 
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Substituting this into fl3.3p gives 

fq{s)+w 



p{t - s, X, q{t) + z)dl.u{x, s) dx 



d^p{t - s, q{s) + w, y)dlu{q{s) + w, s) dy 

q(t)+z 



I rq(t)+z 

- - p{t- s, q{s) + w, y)d'^^u{q{s) + w, s) dy 

+ -/ / p{t- s,x,y)dl-^\{x,s)dxdy. 

^ J —OO J —OO 



We now use the fact that dxp(t, x, y) = —dyp(t, x, y) to rewrite this as 

/q{s)+w 
p{t — s, x, q{t) + z)dl.u{x, s) dx 
'OO 

- ^Pit - q{s) + w, q{t) + z)dlu{q{s) + w, s) 
- l^d'x'^'^i.fli.s) + w,s) / p{t- s,q{s) + w,y)dy 

^ J -OO 

i'q{t)+z i>q{s)+w 



- / p{t- s,x,y)di-^Mx,s)dxdy. (3.4) 

^ J —OO J —OO 



2 

Similarly, for the other partial derivative, we have 
rq{t)+z 

dsC = q'{s) / p{t-s,q{s)+w,y)dl.u{q{s) + w,s)dy 

J —OO 

q(s)+w i'q{t)+z 

dtp{t — s, X, y)dl.u{x, s) dy dx 

i'q{s)+w fq{t)+z 



/q{s)+w rq^tj+z 
/ p{t — s, x,y)dtdl.u{x, s) dy dx. 
■OO J — OO 



By fl3.3p . this becomes 

rq{t)+z 



/q{t)+z 
p{t - s, q{s) + w, y)dlu{q{s) + w, s) dy 
■OO 



l'q{s)+w 

dtC + q'(t) I p{t - s,x,q(t) + z)dl.u{x,s)dx 

-OO 

''q{s)+w i>q{t)^\-z 

p{t — s,x,y)dtdl.u(x,s) dydx. (3.5) 
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Substituting (13. 4 p into the above, and using dtu = {l/2)dlu gives 

^^C = ^P{t - s, q{s) + w, q{t) + z)dlu{q{s) + w, s) 

+ U'(s)9Xg(s) + +w,s)J j p{t - s, q{s) + w,y) dy. (3.6) 

Now observe that, taking i = 0, p(s, t) = ({s, t, 0, 0) — a^. Hence, by (13. 4 p and (II. 2p . 

/<?(«) ]^ 
p(t - s, X, g(t))M(x, s) rfx - -p{t - s, g(s), q{t))u{q{s), s) 
-oo 

+ u{q{s),s)q'{s) p{t - s,q{s),y) dy + - / p{t - s,x,y)dlu{x, s) dx dy, 



oo — oo 



and by ([3lD and ([O]), 



dsp{s, t) = -pit - s, g(s), q{t))u{q{s), s). 



(3.7) 



Differentiating (13. 7p with respect to t gives 



dstpis,t) = -dtp{t - s,q{s),q{t))u{q{s),s) + -9j,p(t - s, q{s), q{t))u{q{s), s)q'{t) . (3.8) 

Part (i) now follows easily from (13. 7p . using the fact that u{q{s), s) is continuous and strictly 
positive on [0, oo), and \q(t) — q{s)\ < C\t — s\ for all s,t E [0, T]. 
Part (ii) will follow from (i), once we show that 



\d,p{s,t) + dtp{s,t)\<C. 



(3.9) 



By dSS]), 



ri{t) 

dsp{s,t) + dtp{s,t) = q'{s) / p{t- s,q{s),y)u{q{s),s)dy 

J — oo 



p{t — s, X, q{t))u{x, s) dx + / / p{t — s, x,y)dtu{x, s) dy dx. (3.10) 

oo J —oo J —oo 



Note that 



q (sj 



q{t) 



p{t - s,q{s),y)u{q{s),s) dy 



\q'{s)\u{q{s),s)P'^^'\B{t - s) < q{t)) < a 

(3.11) 



Next, by the Markov property, u{y, t) = f^p(t — s, x, y)u{x, s) dx. Hence, 



q{s) 



/ p{t — s,x,q{t))u{x, s) dx 



<\q'{t)\u{q{t),t)<C. 



(3.12) 
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Finally, note that dtu{x,t) = (1/2)(9^m(x, t) = (l/2)E^[/"(5(t))]. Therefore, 



OO J — OO 



p{t — s, X, y)dtu{x, s) dy dx 

p{t — s, X, y)p{s, X, z)f"{z) dz dy dx 

= \ f P\B{s) < q[s),B{t) < q{t))f"iz)dz, 



2 



OO J — OO 



which implies 

pit — s, X, y)dtu{x, s) dy dx 



OO •/ — OO 



< ^iiriii- (3.13) 



Applying fl3JT]) . fl3:T2|) . and f lXTSj) to f lXTOj) establishes (jM]), and completes the proof of 
(ii). 

Part (iii) will follow from (13. Sp . once we establish that for |t — s| sufficiently small, 

-C2\t - s\-^/^ < dtp{t - s,q{s),q{t)) < -Ci\t - s\-^/\ and (3.14) 

\dyp{t - s, q{s), qm <C2\t- s\~'/\ (3.15) 

First, note that 

d,p{t,x,y) = --^e-(^-y^'/'^(l ^^"^^'^ ^ 



2v^ V t J t3/2 



Since \q(t)—q{s) \ < C\t — s\, it follows that for \t—s\ sufficiently small, dtp{t— s, q{s) , q{t)) < 
with \dtp{t — s, q{s), q{t))\ < C\t — s\~^^'^, which establishes the lower bound in (I3.14p . For the 
upper bound, note that e"^^^*)"^^''))^/^^*"'') > e"'^l*~''l > e"'^'^. Hence, \dtp(t — s,q{s),q{t))\ > 
C\t — when |t — s| is small enough. 

Finally, we observe that 

\dyP{t,x,y)\ = A=^^^-'^^-'^'"' < \^-y\t-'^^ 



so that (13.150 follows as above from \q(t) — q{s)\ < C\t — s\. This completes the proof of the 
theorem when p is everywhere replaced by p. We now prove the theorem as stated. 

By dO]), p{s,t) = 9{s)e{t)p{s,t), where e{t) = {u{q{t),t))-\ Note that 6 G C^[0,oo) 
and 

dsp{s, t) = e'{s)e{t)p{s, t) + e{s)e{t)d,p{s, t), 
dtp{s, t) = e{s)e'{t)p{s, t) + e{s)e{t)dtp{s, t), 

dlpis, t) = e'{s)e'{t)p{s, t) + e'{s)e{t)dtp{s, t) + e{s)e\t)dsp{s. t) + e{s)e{t)dlp{s, t), 

from which (i), (ii), and (iii) follow. □ 
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Corollary 3.2. Fix T > 0. There exist positive constants Ci,C2 such that 

Ci\t - < E\F{t) - F{s)\^ < C2\t - (3.16) 

for all s,t e [0,T]. 

Proof. Since H- E\F(t) — is continuous and strictly positive on {s ^ t}, it will 

suffice to show there exists Sq > such that the theorem holds whenever |t — s| < Sq. Also, 
in the notation of the proof of Theorem 13.11 since F{t) = 6{t)F{t), where 6 G C^[0,oo) is 
strictly positive, it will suffice to prove the theorem for F. For this, observe that 

E\F{t) - F{s)f = p{t, t) + p(s, s) - 2p{s, t) = 2{p{t, t) - p{s, t)), (3.17) 

where p is given by (13. ip . Hence, 

t 



E\F{t) - F{s)\'^ = 2 I dsp{u,t)du. 

J s 



Applying Theorem I3.1( i) and the fact that J^* \t — u\ du = 2\t — finishes the proof. 
□ 

Corollary 3.3. The process F has a modification which is locally Holder continuous on 
[0, oo) with exponent 7 for any 7 < 1/4. 

Proof. By the Kolmogorov-Centsov theorem (see, for example, [10], Thm 2.2.8), if, for each 
T > 0, there exist positive constants (3,r,C such that 



E\F{t)-F{s)f <C\t-s 



l+r 



for all s,t E [0,T], then F has a continuous modification which is locally Holder continuous 
with exponent 7 for all •y < r/(3. Since F is Gaussian, (13.161) implies E\F(t) — F{s)\p < 
Cp\t-sf\ We may therefore take P = p and r = p/4 — 1. Letting p — j- 00 completes the 
proof. □ 

Corollary 3.4. For each T > 0, there exist positive constants 6,Ci,C2 such that 

-C2\t - s\~^/'^At < E[F{s){F{t + At) - F{t))] < -Ci\t + At - s\~^/'^At, 

for all < s < t < T and At > such that \t — s\ < 6 and At < 6. 

Proof. We write 

/t+At 
dtpis,u) du. 

Applying Theorem I3.1( ii) and the fact that 

-t+At 

It + At - sr^/^At < / |n-sr^/^rfn < It-sr^/^At, 



finishes the proof. □ 
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Corollary 3.5. For each T > 0, there exist positive constants S,Ci,C2 such that 



-C2\t- s\~'^/^AsAt < E[{F{s) - F{s - As)){F{t + At) - F{t))] 

< -Ci\t + At- (s - As)|"3/2^sAt, 

for all < s < t < T, As > 0, and At > such that \t - s\ < 5, As < 5, and At < 5. 
Proof. We write 



E[{F{s) - F{s - As))iF{t + At) - F(t))] = / / dl 

Js-As Jt 



p{u, v) dv du. 



Applying Theorem 13. l( iii) and the fact that 



It + At - (s - As)r^/2AsAt < 



s-As Jt 



t+At 



\V — U\ 



-^'^dvdu < It - sr^/^AsAt, 



finishes the proof. □ 

Remark 3.6. Note that in Corollaries 13.41 and 13. 5[ since dtp{s,t) and d'^^p{s,t) are continuous 
away from {s = t}, we have 

\E[F{s){F{t + At) - F{t))] \ < C\t- s\'^/^At, and 
\E[{F{s) - F{s - As)){F{t + At) - F(t))]| < C\t - sI'^/'AsAt, 

for all < s < t < T, As > 0, and At > 0. Also note that Theorem 13.11 and Corollaries 
[3:211331 all apply to F(t) = u{q{t),t)F{t) as well. 

As an application of these estimates on the local covariance structure of the increments 
of F, we now show that F has a finite, nonzero, deterministic quartic variation. 



Theorem 3.7. // 



Vu{t)= \F{h)-F{t,^^] 



o<t,<t 



where 11 = {0 = tg < ti < t2 < ■ ■ ■ } is a partition of [0, 00) with tj f 00, then for all T > 0, 



lim E 

inKo 

where = supj \tj — tj_i| 



sup 

0<t<T 



Vn{t)-- f\u{q{s),s)r'ds 

Jo 



Proof. We again adopt the notation of the proof of Theorem 13.11 Let 6{t) = {u{q{t),t)) ^ 
and recall that 9 G C^[0,oo) is strictly positive, and that F{t) = 9(t)F(t). Since Vn is 
monotone, it suffices by Dini's theorem to show that 

Vn{t)^- fe{sfds, 
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in for each fixed t > 0. For any stochastic process X, we adopt the notation 
AXj = X{tj) — X(tj^i). We will also write Atj = tj — tj-i- Hence, 

AF,=e(t,)AF,+F(Vi)A%, 
which implies AF^ = 6{tj)^AF^ + Rj, where 



i-i 



i=l 



Thus, 



E 



J2 <cE^ E \F{t,-i)i/^AFf~^ 



0<tj<t i=l 0<tj<t 

Note that Corollary 13.21 also holds for F. Thus, using Holder's inequality, we have 

2 4 



E 



<tj<t i=l '-\o<tj<t / ^0<tj<t ^ 

= ^E( E ^^?)( E Atf-)i^[|F(t,_OP^AF,^(-^)]) 

E ^^?)( E ^tf-'\E\nt,.,)ni\EAF^ 
<c±(j: Ati)( A(f-(w, 



4(4-i)Nl/2 



^4(4-i)Nl/2 



As in the proof of Corollary I3.3[ since F is Gaussian, (KW) implies EAF^^'"'^ < dAt^-'. 
Thus, 

0<tj<t i=l ^0<tj<t ' ^0<tj<t ^ 

as I HI — 0. It therefore suffices to show that 



V Q{t.;)'AF^^- teisfds, 

— TT /n 



0<ti<t 



in as |H| -> 0. 
By dSlZD, 



EAF/ = 2(p(t„ t,) - p{tj_i,tj)) = 2 / dsp{tj - e, t,) de. 



At, 



(3.18) 
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By m, 



dsp{tj - e,tj) = -p{e,q{tj - e),q{tj))u{q{tj - e),tj - e) 



1 



2V27r 



Using |1 - e-^l < x for all x > 0, and \q{t) - q{s)\ < C\t - s| for all < s < t < T, and 
u{q{-), ■) G C^[0, oo), this gives 



dsp{tj - e, tj 



Substituting this into fl3.18p gives 



1 



-1/2 



2V27r 



u{q{tj),tj] 



EAFf-J^Atfu{q{t,),ti, 



< CAtf\ 



Since EAF/ = ?,{EAFff and e{t) = {u{q{t) , t))-\ this implies 



e{tjfEAFf - -AtjOitj) 



EAFf + \I^Atf9{t,r' 



< CAt'j. 



Note that 



0<tj <t 

It will therefore suffice to show that 



y -AtMtjf - f e{sfds. 



E 



J2 m,Y^Ff-e{t,YEAFf) 



0<ti<t 



^0, 



as |n| 0. 

For this, we write 



E 



J2 mfAFf - eit.fEAFf) = Y,0{Ufe{t,)\E[AFtAFf] - {EAFt){EAFf)). 



0<tj<t ij 

If X and Y are jointly normal with mean zero and variances aj^ and ay, then 

E[X'Y^] = 44(9 + 72ply + 2Apy), 
where px,Y = {o-xO'y)~^E[XY]. Hence, 

lEiX^Y""] - {EX^){EY^)\ < Calal-\E[XY]\^. 
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Therefore, 



E 



0<tj<t i,j 



By Holder's inequality, \EAFiAFj\''^ < C At]''^ At]'"^ . Thus, it will suffice to show that 

Aty^At'/^\EAFiAFj\^ 0, 

tj-ti>2\U\ 

as |n| 0. 

If tj —ti> 2|n|, then tj-i > ti. Hence, by Corollary 13.51 and Remark 13.61 

\EAFiAFj\'^ < CAt^At^lt^.i -t.|-3, 
where C depends only on T. Since \tj-i — tj| > |H| > Atk for all k, this implies 

\EAFiAF,\^ < CAtJ/'Atf - U\-^/' < C\Uf/'Aty'Atf\t,.i - U\-^/\ 

Hence, 

J2 Atj/'Atf '|EAFiAF,f < C|H|3/4 ^ At^At, - t^r^/^ ^ 0, 

tj-ti>2\U\ tj~ti>2\U\ 

since \u — v\~'^l'^ du dv < oo. □ 

Although the local properties of F are very similar to 5^/^, the global properties need not 
be. One simple observation along these lines is that E\B^{t)\'^ = t"^^ , whereas E\F{t)\'^ > Ct, 
for some constant C that depends on the initial distribution /. Indeed, by (11. 3p . we have 
E\F{t)\^ = {a- a^)\u{q{t),t)\-^, and 



uix,t) = -= / fiy)e"^''-y> "'dy< 



2711 J ^/2T\t 

In other words, for t large, F{t) has a variance which is comparable to Brownian motion, 
rather than fBm. 

We will illustrate some other global properties of F with a particular example of an initial 
distribution. As noted in Section [H if -B(O) ~ -/V(0, 1) and j(n) = \ {n + 1)/2J, so that Qn 
is the median, a = 1/2, and q = 0, then the result in [18] implies that Fn{-) =r- X{- + 1) in 
C[0, oo), where the covariance function of X is given by (11. 4p . We will illustrate the global 
properties of F in this case by deriving some of the global properties of X. 

First, note that fBm has a self-similarity scaling property. Namely, (ct) and B^{t) 

have the same law, as processes. Based on this, we might expect X to have a similar scaling 
property with exponent 1/4. However, X obeys the Brownian scaling law with exponent 
1/2. This can be seen immediately by direct inspection of (II. 4p . 

Second, note that fBm is a long-memory process. In particular, 

rnin) = £;[fi^(l)(5^(n + 1) - 
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r 



decays only polynomially. It is well-known that r//(n) ~ 2H{2H — l)n'^^ ^, where a„ ~ 6„ 
means an/bn — > 1 as n — )■ oo. Based on this, we might expect 

r{n) = E[X{l){X{n + l) - X{n))] (3.19) 

to satisfy r{n) ~ some constant C. However, this is not the case. 

Proposition 3.8. If r{n) is given by f l3.19p . then r{n) ~ — (l/6)n~^. 
Proof. The proposition follows easily by observing that 

-(n) = y/n + 1 sin~^ ( —=] — \pn sin~^ f ) 

and using the Taylor expansion sin~^ x = x + + 3x^/40 + 0(x''). □ 

We see then that X is a process which behaves locally like fBm with H = 1/4, but whose 
increments are more weakly correlated than those of fBm. Another example of such a process 
is sub-fBm (see Bojdecki et al [3], for example). For sub-fBm with H = 1/4, the increments 
decay at the rate n~^l'^ . The decay rate of n"^ in Proposition 13.81 has only been established 
here for the case when the initial distribution of the particles is a standard Gaussian, and an 
investigation of the relationship between the initial distribution and this exponent is beyond 
the scope of the current paper. However, we may presently note that F is not a sub-fBm 
for any initial distribution /, since the variance of a sub-fBm is always proportional to t^^, 
whereas E\F{f)^ > Ct, as we observed above. 

4 Outline of proof of tightness 

Our primary tool for proving tightness is the following version of the Kolmogorov-Prohorov 
tightness criterion. 

Theorem 4.1. If {Zn} is a sequence of continuous stochastic processes such that 
(i) sup„>„g P{\Zn{fS)\ > A) — 7- as A — oo, and 

(li) sup„>„^,P(|Z„(t) - Zn{s)\ > e) < Cre-'^lt - s|^+^ whenever < e < 1, T > 0, and 
0<s<t<T, 

for some positive constants no, a, P, and Ct (depending on T), then {Zn} is relatively 
compact in C[0, oo). 

When we apply the above theorem to our processes {Fn}, Condition (i) will automatically 
be satisfied due to the already established convergence of the fdds. Our main theorem will 
therefore be proved once we establish that {Fn} satisfies Condition (ii) of Theorem 14.11 The 
remainder of the paper will be devoted to proving this. For this, we will need to study the 
distribution of the difference of two quantiles. Unfortunately, the difference of the quantiles is 
not the quantile of the differences. The stark non-linearity of the quantile mapping creates 
substantial technical difficulties. In the next section, we begin by using conditioning to 
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connect the probabihty being estimated to a certain random walk. With an eye to future 
apphcations, we construct this connection in abstract generahty. In Section [6l we will return 
to the specific sequence {Fn}. 

Before moving on to the details of the proof of tightness, we pause here to present a 
general overview of the key steps to be taken in the remainder of the paper. Let Bj = Bj — q, 
B = Bi, and = Qn — 1 = Bj(^n):n- To verify Condition (ii) of Theorem 14. ![ we will aim 
to show that 

P(|F„(t) - >e)< Cpe-P\t - s\p/\ (4.1) 

for any p > 2. We begin by writing 

P{\Fn{t) - F^{s)\ >e) = P{mt)-Qnis)\ > n-'/'e). 

Heuristically, the order of magnitude of \Qnit) — should be less than that of 

\B(t) — B{s)\, which is |t — Consequently, it is easy to estimate this probability when 
n~^/^£: > |t — for some A > 0. (This is what we call the "large gap regime".) This 

routine estimate is carried out in the appendix, and (14. ip is easily verified in the large gap 
regime. 

We next turn our attention to the "small gap regime", when n~^/^£ ^ \t — for 
some A > 0. We begin with the fact that 

Pmt)-Q^s) < -n^l\) = E\P{QSt)-ms) < -n'l^e \ QM)\- (4-2) 

(The estimates for P(Q^(t) — Q^is) > n^^^'^e) are nearly identical.) To deal with the right- 
hand side of (14. 2p . we apply the results from Section |5l wherein we study, in generality, 
the conditional distribution of the difference of two quantiles. There we find that this 
conditional probability can be well-approximated by certain probabilities connected to sums 
of iid random variables, and that sufficient estimates on these probabilities can be derived. 
When we apply Theorem 15.11 to the problem at hand, we find that if 

qi{x, y) = P{B{t) > q{t) + x + y\ B{s) < q{s) + x), 
q2ix, y) = P{B{t) < q{t) +x + y\ B{s) > q{s) + x), 

and 

^f-.ni^^y) =P[ $^lW<gi} < l{C/,<<?2} 1) 

^ 1=1 i=j+l ^ 

where {Ui\ are iid and uniform on (0, 1), then 

P{Qnit)-Qn{s) < -n'l\ I g„(s)) < y,|(„^^jg„(.),-n-V2^). 

By (14. 2 p and the fact that ipj.^{x, y) < I, this gives, for any A' > 0, 

P(Qnit)-Qr.is) < -n^l\) < sup cpf {x,-n-'^'e) + Pi\Q,Xs)\ > K). 

\x\<K ^ ' 

Since F„ has well-behaved tail probabilities (see Proposition IA.2I and Remark lA.Sp . we 
may choose K so that P(|(5„(s)| > K) < Cpe"^\t — s\^^^. (In fact, the correct choice is 
K = n~^^'^e\t — s\~^^^.) The problem thus reduces to showing that 

sup ^7h.„(x, -n-'^h) < C,e-^\t - sf\ (4.3) 

|x-|<X 
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For this, we again appeal to the general results in Section [51 This time applying Theorem 
15.21 we find that 

where 

cr(x, y) = aqi{x, y) + {1 - a)q2{x, y), 
y) = aqi{x, y) - {1 - a)q2{x, y). 

To show that this leads to (14. 3p . we use Taylor expansions for a and /i which are developed 
in Proposition 16.11 and Corollary 16.21 and which rely on some of the computations in the 
proof of Theorem 13. 11 When all of this analysis is finally complete, we will have established 
(14. Sp . and therefore (14. ip . in the small gap regime. 

Lastly, we must deal with the "medium gap regime", in which \t — sl^/^"*"^ < n~^/'^e < 
\t — s|^/^~^ for some A > 0. This case is treated by making minor modifications to the proof 
for the small gap regime. In making these modifications, however, the result is weakened, 
and we are only able to prove that, for parameter values in the medium gap regime, 

P(|i^„(t) - >e)< C{e-'\t - s\'/'-'^y. (4.4) 

Although we conjecture that the sharper bound (14. ip does in fact hold in all regimes, the 
weaker bound (14. 4p which we are able to prove is still sufficient to establish tightness. 



5 A random walk representation and estimate 

Let {Un}'^=i be an iid sequence of random variables, uniformly distributed on (0, 1), and for 
'"i, r2 G (0, 1), define 

^J^„(ri,r2) = Pf ^l{C/,<ri} < ^{U^<r2}j, 
^ i=l j=j + l ^ 

i^f:.ni^Ur2) = Pf ^l{C/,<r-i} < ^{U^<r2}]■ 

Also, let tp>^ = 1 - i)f.^^ and = 1 - V^<„. 

Let X and Y be real-valued random variables such that (x, y) t— ?■ P{X < x,Y < y) is 
continuous. Define 

gi(x, y) = P(Y > X + y \ X < x), 
q2{x, y) = P{Y < x + y \ X > x), 

and let (p*.j^{x, y) = ip*.j^{qi{x, y), q2{x, y)), where -k may be any one of the symbols <, <, >, 
or >. Note that ^f,^{x,y) = V'(";-j+i):„(92, gi) and ^f,^{x,y) = ^f^_j^^yjq2,qi)- 
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Theorem 5.1. // Yn)} is a sequence of iid copies of (X, Y), then for all ?/ G M, 

iffAXr.n. y) < P{Yr.n - < y \ X,,„) < (/.|„(X,:„, y), (5.1) 
almost surely. Consequently, 

<^>„(X,,„, y) < PiYj.,^ - > y \ X,,„) < v%{X,..n, y), (5.2) 

almost surely. 

Theorem 5.1 estabhshes a connection between the conditional distribution of the 
difference of two quantiles and probabilities connected to sums of iid random variables. To 
give some intuitive sense to this connection, let us consider the following heuristic derivation 

of (EU). 

We are interested in estimating P(Y^:„— Xj;„, < y \ Xj,n = x) = PiYj.n < x+y \ Xj.n = x). 
Let us consider Xi, . . . ,X„ as representing the locations on the real line of some particles 
at time s, and Yi, . . . ,Yn as representing the locations of those same particles at some later 
time t > s. We are given that the j'-th leftmost particle at time s is located at position x. 
Conditioned on this information, we know the following. At time s, there is one particle 
located at x, there are j — 1 iid particles located in (— oo, x), and there are n — j iid particles 
located in (x, oo). The event {Yj,n < x + y} will occur if and only if at least j particles end 
up in (— oo, x + y) at time t. 

Each particle which is in (— oo,x) at time s has probability 

Pi = 1 - gi = P{Y < x + y \ X < x) 

of ending up in the target interval (— oo,x + y) at time t. Therefore, we may represent the 
number of particles from (— oo,x) which end up in the target interval by ^^=1 l{c/i<pl}• 
Similarly, the number of particles from (x, oo) which end up in the target interval is 
represented by YlLj+i l{c/,<g2}- 

We are therefore interested in computing the probability that 

j-l n 

1=1 i=j + l 

j~l n 

i=l i=j+l 
j-1 n 

1=1 i=j+l 

which happens if and only if Yl^i=i ^{u^<ql} < Yl^=j+i ^{Ui<q2}- This probability is exactly 
ipf.^{qi,q2) = (pf.n{x,y). In fact, the true probability is even larger than this, since the 
particle at x may itself end up in the target interval (— oo, x + y) at time t. If that happens, 
then we only need YllZl l{c/i<gi} < Zir=i+i '^{Ui<q2}y and the probability of this is (ff.^{x,y). 
Hence, the conditional probability of interest is sandwiched between ^9*:^ and (pj.^. 
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Proof of Theorem 15.11 By taking complements, ( 15. ip and (15. 2p are equivalent. We will 
prove (15. 2p . For this, we will first establish that for any a, 6, y G M with a < b, 



E[^f:ni^j:n,y)UXr.ne(a,b)}] ^ ^i^Pn " ^r-n > V.Xj.n E (a, 6)) 

< ^bL(^j:n,y)l{X,:„e(a,6)}]- (5-3) 



Let us begin with the upper bound in (15. 3p . 

First, by taking ai ^ a and bi b if necessary, we may assume that P{X < a) > 
and P{X > b) > 0, so that g^, ^ff-n, and where ^{x) = P{X < x), are all well-defined, 
continuous, and bounded on {{x,y) G [a,b] x M}. Now fix £ > and x G [a,b — e] and let 
N = #{i: XiE {x, X + e)}. Then 

P(y,:„ - X,,„ > y, G (x, X + 5)) < P(iV > 2) 

+ P(r,,„ > X + #{z : X, < x} = J - 1, #{z : X, > X + e} = n - J, AT = 1). 

Thus, 

P{Yj.,n - Xj,n > y, Xj,n G (x, X + e)) 

- (2) ^ + + ^' (j) ^^^^'^^ > ^ + ^' (5-4) 

where 

A = {Xi < X for all i < j} (1 {Xi > x + e for alH > j} fl {Xj G (x, x + e)}. 

Note that 

P(r,;„ > X + y, A) 

= Pmt:Y,<x + y}<j,A) 
<Pmt^j:Y,<x + y}<j,A) 
i-i i-i 

= ^ J]P(#{z<j <x + y} = £,#{z>j <x + 2/} = m-£,A). 

i=0 m.=l 

By independence, this gives 



P{Jy,n > X + A) 

< 

^=Q m=£ 



JP(>^<x + y,X<x)^P(F>x + i/,X<x)^-^-^ 



m — i 



^ 1 P(F < X + X > X + e)'"-^P(F >x + y,X>x + 



P(X G (x,x + 5)). 



21 



If we define qi{x,y) = qi{x + e,y — e), pi = 1 — qi, pi = 1 — qi, and $ = 1 — $, tlien tliis 
becomes 

i-i j-i 

j — l—£~n—j—m+£~m—i 



■ ^{xy-^^{x + eY-^P{X G (x, X + e)). (5.5) 



Compare tliis witli 

i-i 



^ i=l i=j+l ' 

£=0 m=i ^ ^ ^ ^ 



Finally, partition (a, 6) into subintervals of size e and apply (15.41) and (15. 5p . Let e — >■ and 
apply dominated convergence. By (15. 6p . this gives 

P(F,,„-X,,„, > G (a,6)) < ^V|;n(a;,i/)'^(a;r'^(a^)"^' ^^*(a;)- 

By Lemma [A. 41 this completes the proof of the upper bound in (15. 3p . 
The lower bound in (15. 3 p is proved similarly. As before, 

P(Y,-;„ - Xj.,n > y, Xj,n G (x, X + e)) 

> P{Yj.,n > X + e + #{z : < x} = J - 1, #{z : X, > X + e} = n - J, AT = 1), 

which gives 

P{Yj.,n - Xj.,r, > y, Xj.,n G (x, X + s)) > J Q P(Yj.,r, >x + e + y,A) 

n 
J 



> J ( ) y^j:Y,<x + e + y}<j-l,A). 



Note that 



P{#{i ^j:Y,<x + e + y}<j-l,A) 

= ^ ^ P(#{2 <j:Y,<x + e + y} = i,i^{t>j:Yi,<x + e + y}=m-i,A). 



1=0 m=i 
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If we now define qi{x, y) = qi{x, y + e), q^i^x, y) = qi{x + e,y), % = I - and Pi = 1 - 
tfien tfiis gives 

P{Yj,n - Xj.,n > y, Xj,n G (x, X + e)) 

\ i-2 i-2 / . ^ ^ / 

sr^ / 



■ <l>(x)^"i$(x + £)"--'P(X G (x, X + e)) 

We compare tliis witli 



^ i=l i=i+l ''^ 

i-2 i-2 . i-1 n N 

^=0 m=£ ^ 1=1 i=i+l ^ 



e=0 m.=i 



and the remainder of the proof is the same. 

Using (15.31) . we now prove the upper bound in (15.21) . Let 

^ = P(Yj.,n - Xj.,n > y I Xj.,n) - ifJ-A^j-.n, y) , 

SO that > 0} = {Xj.,n e B} for some Borel subset B cR. Fix e > 0. There exists Bq cR 
such that Bq is a finite, disjoint union of open intervals, and 

\P{X^.,^eB)-P{X,.,r,eBo)\<e. 

(See Proposition 1.20 in [6], for example.) Hence, by (15. 3p . 

= P{Yj.,n - Xj.,n > y, Xj.,n G B) 

< P(Yj.,n - Xj.,n > y, Xj.,n eBo)+e 

< E['Pf:niXr-n,y)'^{x,..,,eBo}] + £ 
<E[^f^^{X,.,r^,y)l{x,..„eB}]+2e. 

Therefore, i?[^l{^>o}] < 2e. Letting e -> shows that ^ < a.s., completing the proof. The 
lower bound in (15.21) is proved similarly. □ 



With Theorem 15.11 we have accomplished the first step of connecting our quantiles to a 
random walk. The second step, given by the next theorem, is to derive an estimate for this 
walk. 

Theorem 5.2. Fix a G (0, 1). Let 

a = a{ri, r2) = ari + (1 - a)r2, 
/i = yu(ri, rs) = ari - (1 - a)r2. 
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and suppose j{n)/n = a + o{n '^^'^). Then for each t > 1, there exist constants C > and 
riQ G N such that, for all n > no, 

V^j-(n):n(n,r2) < C 



whenever fi > 0. Note that C does not depend on ri or r2. 
Proof. First note that if nari < 2, then 



Hence, we may assume that nari > 2. 

Now, since j{n)/n = a + n~^^'^an, where a„ — )■ 0, we may write 

V',-(n):Jn,r2) = - ^u) + bn< -/i) < P{\n~\^L - ^u) + bn\ > /i), 

where 



i=l 
n 



i=j{n)+l 

bn = n~^/^a„(ri + - ra~Vi. 

By Chebyshev's inequahty. 
Note that 

^ -1/2 



a A (1 — a) ' 

which imphes |fenP^ < Cn^'^a'^'^ . It will therefore suffice to show that 

E\^l\'^ + E\^u\'^ <Cinay. 
By Lemma [A.51 there exists and C such that n > uq implies 

E\^L\'-<C{U{n)r^yy{J{n)r^)). 

By making uq larger if necessary, we may assume na/2 < j{n) < ?)na/2. Since nari > 2, 
this gives j{n)ri > 1, so that 

E\^l\^^ < C{j{n)riy < CinariY < C{nay. 

Similarly, for n sufficiently large, 

E\^u\'^ < C{{{n - j{n))r2r V {{n - j{n))r2)) 
<C((n(l-a)r2rvl) 
< C((n(l - a)r2y V (nari)") < C{nay, 
completing the proof. □ 
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6 Parameter estimates and gap regimes 

Let us now return to the specific assumptions of our model, as outlined in the beginning of 
Section O Fix < s < t. We shall adopt the notation and definitions of Section |5l in the 
special case that X = B{s) — q{s) and Y = B{t) — q{t). In this case, 

qi{x, y) = P{B{t) > q{t) + x + y\ B{s) < q{s) + x), 
q2{x, y) = P{B{t) < q{t) + x + y\ B{s) > q{s) + x). 

Let us also define 

^(x, y) = P{B{t) > q{t) +x + y, B{s) < q{s) + x). (6.2) 

Our objective is to verify Condition (ii) of Theorem 14.11 Ideally, we would like to show that 

P{\F^{t) - >e)< C,e-^\t - sY'\ (6.3) 

In the end, we will actually show something slightly weaker, although our final estimate will 
be sufficient to verify the conditions of Theorem 14.11 Our approach will begin by conditioning 
on Fn{s), so that we may apply Theorem 15.11 We will then use Theorem 15.21 to obtain the 
specific bound we need. Implementing this strategy will require precise estimates on the 
function qi and g2, in terms of x, s, and t. These estimates will come from Taylor 
expansions. We therefore begin with a Taylor expansion for ^ . 

Proposition 6.1. Fix T,K > 0. There exists a constant C such that for all < s < t < T 

and all x, y satisfying \x\ + \y\ < K , 

^{x, y) = *(0, 0) - ^u{q{s), s)y + ^-^u{q{s), s)y^ + R, (6.4) 

where 

\R\ < C {{\x\ + \y\){6'" + \y\ + r^^j^p) + 6-'/^') , (6.5) 

and 6 = t — s. 

Proof. Fix T,K > 0, s,t e [0,T], and x, y G M such that s < t and |x| + \y\ < K. Let 
5 = t — s. In what follows, C shall denote a constant that depends only on T and and 
may change value from line to line. 
By Taylor's theorem, we may write 

^(x, y) = *(0, 0) + a,.^(0, 0)x + a,^(0, 0)y 

+ ^92^(0, 0)x2 + a,9,*(0, 0)xy + ^9J^(0, 0)1/2 

+ \dl^!{ex, ey)x' + Idldy^iOx, ey)x^y + \d,dl^>{ex, ey)xy^ + ^aj^(^x, ey)y\ (6.6) 
2 2 6 
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where 9 e (0, 1). Let = (27r)~i/2 j^^ g-j/Vs dy a.nd ^ = I - $. We first establish that 
for all integers i > and j > 1 , 



q{s)+x 



X + ?/ + - z 



9>(;2, s) dz, (6.7) 
9Mg(s) + x,s) + 9:+ia^-ivl;. (6.8) 



If i = 0, then (16.71) follows directly from the definition of \E', (16. 2p . Differentiating (16. 7p gives 



.u{q{s) + X, s) 



<j(s)+x 



X + 1/ + ^(t) - Z 



dl.u{z, s) dz. 



Applying integration by parts shows that (16.71) holds for i + 1. By induction, this proves 
(16. 7p . For (16. 8p . let j = 1 and let z > be arbitrary. By (16. 7p . we have 



_ f X + y + q{t) 

-00 

q(s)+x 



5V2 



- / x + y + q{t) - z 
' 51/2 



y + q{t) -q{s) 
51/2 



^ ( y + q{t) - q{s) 
' 51/2 



9*u(z, s) dz 

dl.u{z, s) dz 

d'xu{q{s) + X, s) 
diu{q{s)+x,s)+di+^^ 



which is (16. Sp . Now assume (16. 8p is valid for some J > 1 and all i > 0. Applying dy to both 
sides of (16. 8 p immediately shows that (16. 8 p is valid for j + I and any i > 0. By induction, 
this proves (16. 8p . 

Now, by (16. 7p . we may write, for any i > 1, 



di^{x,y) 



q(s)+x 



dl.u{z, s) dz 



dl u{q{s) + x,s) - C{s,t,x,x + y), 



q{s)+x i>q(t)+x+y 



pit — s, z, w)dl.u{z, s) dw dz 



where ( is given by (13. 2p . We shall adopt the convention that ^u{x,t) := P{B{t) < x), so 
that the above equality is, in fact, valid for all i > 0. Hence, 



dl.'^{x,y) = dl u{q{s) + x,s)+ / dsC{r,t, x, x + y) dr - ({t,t, x, x + y). 



Since 



C{s,t,w,z) 



q{s)+w 



P^(E(t -s) < qit) + 2)9>(x, s) dx, 
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it follows that 



C{t,t,w,z) 



q{t)+{whz) 



9Xx, t) dx = d'~^u{q{t) + {wAz),t). 



Therefore, 



y) = d'-^u{q{s) + X, s) - d'-^u{q{t) + x + A 0), s) + / dsC{r, t, x, x + y) dr. 

J s 

By the mean value theorem, since |g(t) — q{s)\ < C\t — s\ = C6, we have 

\dl-\{q{s) + x,s)- dlr\{qit) +x + iyAO),s)\<C{S+\y\). 
By fl3.6p . \dsC{r, t, x, x + y)\ < C\t — rl^^^"^, and so we obtain 



\dl^{x,y)\<C{6'/' + \y\), 



(6.9) 



for any i > 0. (Here we have used the fact that since 6 < T, there exists C such that 
6 < C6^/^ for all 6<T.) 

We next consider the derivatives with respect to y. By (16. 8p . 

'g(t) - q{s) 



^2,^(0,0) = -$ 



u{q{s),s) + 8^^(0,0). 



By O, 19^.^(0,0)1 < C6^/\ Also, \q{t) - q{s)\ < C6 and <l>(x) = l/2 + 0(|x|). Hence, 



dy^{0,0) + -u{q{s),s] 



Similarly, 



and 



|<9.a,^(o,o)| 



g(t) - g(s) 



51/2 



< c, 



92^(0,0) = -r^/^<i> 



-/ / g(t) - q{s] 



u{q{s),s) + d,dy^{0,0). 



Since $ (x) = -(27r)-i/2e 



1/2 -xV2 



5,^(0,0) 



^1/2 / -ViV"/' "7 ' -x-y 

(27r)-V2(i + Odxp)), this implies 
1 



V2^ 



u{q{s),s] 



< C. 



For the third-order partial derivatives, we have 



\dldy^{x,y)\ 
\d,d'y^{x,y)\ 
\d'^{x,y)\ 



51/2 



d'M<l{s) + x,s) + dt^{x,y) 



51/2 



51/2 

<Ci5-'%\ + 6~'/% 



dxu{q{s) + x,s) + dldy^!{x, y) 



u{q{s) + X, s) + dxdl^{x, y) 
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where the last inequahty follows from $ (x) < C\x\. 

Substituting all of this into (16. 6 p gives (16. 4p . where R satisfies 

\R\ < C{5^'^\x\ + 5^'^\y\ + 5^'^\x\^ + \xy\ + \y\^ 
Since X and y are restricted to a compact set, we may simplify this to 

\R\ < C{5'/M + ^'"\y\ + \xy\ + \y? + 5-"'\xy'\ + r^/^^r + 5'"'\y?). 

which is precisely (16.50 . 



□ 



Corollary 6.2. Recall qi,q2 given by (16. ip . FixT,K > 0. There exist constants 6o and C 
such that for all < s < t < T and all x, y satisfying \x\ + \y\ < K and \x\ < 6o, 



aq,{x, y) = ^(0, 0) - s)y + — 1=-^"^-^ -^-^ 



(1 - a)q2{x, y) = ^(0, 0) + Ui{q{s), s)y + — 1=-^'-^"^ -^-^ 



M(g(s),s)?/^ + i?2, 



(6.10) 
(6.11) 



where Ri, i?2 both satisfy (I6.5p . and 5 = t — s. 
Proof. By ([23]) and fl^ . 

\^ix,y)\<Ci5'/' + \y\+5~'%f) + \R\, 

where /2 satisfies dSSD- By (EI]) and flO]) . 

^ ^(x,y) 
1i{x.y) p(^B{s) <q{s) + x)' 

Note that P{B{s) < q{s) + x) = a + r{x, s), where |r(x, s)| < C\x\. Hence 

r(x, s) 



(6.12) 



\aqi{x,y) - ^(x,?/)| 



^(x,y) < | ^(x,y). 
— o |x| 



a + r[x, s 

If X is sufficiently small, then using (I6.12p . we have 

\aq^{x,y) - ^ix,y)\ < CiS'/^\x\ + \xy\ + S-'/^\xy^\) + \R\. 
By (16. 4p . this gives 



aqi{x, y) - ^'(0, 0) + ]^u{q{s), s)y - ^-^==M(g(s), s)y^ 



< C{6^/^\x\ + \xy\ + 6-^/^\xy^\) + 2\R\. 



Since this is bounded above by the right-hand side of (16. 5p . this completes the proof of 
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For (16.1ip . let B{t) = —B(t), and let q(t) be the (1 — a)-quantile of the law of B(t), so 
that q{t) = —q{t). Let u{x,t) be the density of B{t), so that u{x,t) = u{—x,t). Define 



and 



^(x, y) = P{B{t) > q{t) +x + y, B{s) < q{s) + x), 

y) = P{B{t) >q{t)+x + y\ B{s) < q{s) + x) 
= P{B{t) < q{t) -x-y\ B{s) > q{s) - x) 
= q2{-x,-y). 



Hence, by (1630|) . 



(1 - a)q2{x, y) = {1- a)qi{-x, -y) 

= §(0, 0) + ^u{q{s), s)y + -^u{q{s), s)y^ + 

= 5(0, 0) + ^u{q{s), s)y + -^u{q{s), s)y^ + i?2, 

where R2 satisfies (16.51) . To complete the proof, we observe that 

5(0,0) = P(i?(t)<g(t),i?(s)>g(s)) 

= P{B{t) < q{t)) - P{B{t) < q{t), Bis) < g(s)) 
= P{B{s) < q{s)) - P{B{t) < qit), B{s) < q{s)) 
= P{B{t) > q{t),B{s) < q{s)) = ^{0,0), 

which gives (I6.1ip . □ 

We are now ready to establish (16. 3p and complete the proof of our main result. Theorem 
11.41 Recall that F„ = n^^'^{Qn — q)- Hence, the event {|F„(t) — Fn{s)\ > e} is precisely the 
event that the process Qn — q traverses a gap of size n~^^'^€ in the time interval [s, t\. In that 
same time interval, each Brownian particle will traverse a gap of size roug hly \t-s\^/^. We 
therefore divide our proof into three separate cases, or parameter regimes: n~^^'^e ^ |t — s|^/^ 
(the large gap regime), n~^^'^€ <^ \t — s]^^"^ (the small gap regime), and n~^^'^e ~ |t — sl"*^/^ 
(the medium gap regime). 

The proof in the large gap regime (Lemma 16. 3p uses only crude estimates and requires 
none of the machinery we have so far developed. The proof is entirely standard and is given 
in the appendix. The proof in the small gap regime (Lemma 16. 4p requires the most care. We 
follow the strategy outlined earlier, utilizing all of the tools developed previously. The proof 
in the medium gap regime (Lemma 16.50 is a minor modification of the proof of Lemma 16.41 
In making this modification, however, our estimate loses some precision, and the result in 
the medium gap regime is, in fact, slightly weaker than (16.31) . 

Lemma 6.3. Fix T > 0, A G (0, 1/2), and p > 2. There exists a positive constant C such 
that 

P{\F^{t) - F„(s)| >e)< C{s^^\t - s\'/y, 
whenever s,t E [0,T], < e < 1, and n^^/^e > \t - s|^/^~'^. 
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Lemma 6.4. Fix T > 0, A G (0, 1/2), and p > 2. There exists a positive constant C and 
an integer uq such that 

P{\F^{t) - > e) < C{e~^\t - s\''y, 

whenever n > uq, s,t E [0, T], < e < 1, and n~^^'^e < \t — sl^/^"*"^. 

Proof. Fix T > 0, A G (0, 1/2), and p > 2. Without loss of generality, assume s < t and 
define 6 = t — s. Note that it is sufficient to prove that there exists a constant 6^ > such 
that the lemma holds whenever 6 < 6o. 

Let no be as in Theorem 15.21 with t = p/2. Define Bj = Bj — q and = Qn — q. 
Note that Q„ = Bj^ny,n and F„ = n^^^Q^. Let K = n'^/'^ed-^''^ and y = n'^^'^e. Then by 
Theorem I5.H 

P(F„(t) - F^is) < -e) = P(Q^{t) - Q^s) < -y) 

= E[P{Q^{t)-Q.M<-y\Qn{s))] 
<EWlnyMn{s),-y)] 

< sup ^%)Jx, -y) + P(|Q„(s)| > K). (6.13) 

\x\<K ^ ' 

By making larger, if necessary, we may apply Proposition IA.2I and Remark IA.3l to obtain 

P{\Qn{s)\ >K) = P(|F„(s)| > e5~^'^) < ^(er^/^)-^ = C{e~^\t - s\^'^Y. (6.14) 

For the other term, fix x G [—K^ K] . Following the notation of Theorem 15. 2[ define 

a{x, y) = aqi{x, y) + {1 - a)q2{x, y), 
/i(x, y) = aqi{x, y) - {1 - a)q2ix, y). 

By Corollary 16. 2 [ 

Kx, -y) = u{q{s), s)y + R^, 

where P^ satisfies (E3]). Note that \x\ < K = n~^l'^Eb~^l^ and y = vT^I'^e < 5^^^+^. In 
particular, y < K and y < 5^^^. Therefore, (16. 5p simplifies to 

Hence, since s t-> u{q{s),s) is bounded below on [0,T], there exists 5q such that 5 < 5q 
implies /i(x, —y) > Cy. Similarly, 

a(x, -y) = 2^(0, 0) + -^u{q{s), s)y^ + P., 

where \R^\ < C{5^/^ + 5'^^)y < Cy < CS^I'^. By (EH]), this implies \a[x, -y)\ < C6^/\ 
We now apply Theorem 15.21 with t = p/2 to obtain 
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Substituting f l6J[5|) and f lCTl) into fl67[3|) gives 

P(F„(t) - < -e) < Cie-'\t - s\'/y. 
To complete the proof, we follow as in f l6.13p to obtain 

P{F4t) - F4s) > e) < sup ^f^^y^{x,y) + Pm,{s)\>K), 

\x\<K ^ ' 

and apply (I6.14p to the second term. For the first term, recall from Section |5] that 
Since {n — j{n) + l)/n = (1 — a) + o(n~^/^), we may apply Theorem 15.21 to obtain 



(y{x,y) 

2 



p/2 



where 

y) = {l- a)q2{x, y) + aqi{x, y) = a{x, y), 
y) = {l- a)q2{x, y) - aqi{x, y) = y). 

The estimates preceding f l6.15p can now be used, as before, to give 

which completes the proof. □ 

Lemma 6.5. Fix T > 0, A G (0, 1/8), and p > 2. There exists a positive constant C and 
an integer uq such that 

P{\F,{t) - >e)< C{e-'\t - sl'/^'-'^Y, 

whenever n > uq, s,t e [0,T], < £ < 1, and \t - s|^/^+^ < n^^^'^e < \t - 

Proof. Fix T > 0, A G (0, 1/8), and p > 2. Without loss of generality, assume s < t and 
define 6 = t — s. Note that it is sufficient to prove that there exists a constant > such 
that the lemma holds whenever 6 < 6o- 

Let no be as in Theorem 15. 2[ with r = p/2. Let K = n~^^'^e5~^^^ and y = rT^^'^e. Also 
define e = n^/^^Vs+A < ^ and y = n-^l'^e = (5V2+a_ in (EH and flCTl) . 

P(F„(t) - F„(s) < -e) < P{F^{t) - F„(s) < -e) 

< sup Jx,-y) + C{e-'\t-s\'/'r. 

\x\<K ^ ' 

By the estimates following f l6.14p . ^{x, —y) = u(q{s), s)y + R^, where 
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Since A < 1/8 and s u{q{s),s) is bounded below on [0,T], there exists 60 such that 
S < So implies fi{x,—y) > Cy. Similarly, \a(x,—y)\ < Cb^l'^. As in (16.151) . this gives 
^\n)A^^ -y) < C{rr^y^'^b^l'^Yl'^. Note that y = 5^/2+^ > rT^I'^eb'^^. Hence, 

^InyA^^ -V) < C{e-'6'/'-'y/' = C{e-'\t - .jV^-^^f. 

As in the second half of the proof of Lemma 16.41 the bound on P{Fn{t) — -F„(s) > e) is 
proved similarly. □ 

Proof of Theorem 11.41 By Corollary 12.21 it will suffice to show that {F^} is relatively 
compact in C[0, 00). By Theorem 14.11 we need only verify that 



sup P(|F„(t) - F„(s)| >e)< - s 

n>no 



1+/3 



Taking A = 1/16 and p = 9 in Lemmas 16.31 16.41 and 16.51 shows that 

P{\F4t)-F4s)\>e)<Ce~'\t-sf/', 
for all n > riQ, which completes the proof. □ 
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A Appendix 

In this section, we have collected together the more routine and straightforward proofs, so 
as not to distract from the essential techniques of the main paper. We begin with the proof 
of Lemma II. 1[ 

Proof of Lemma II. IL Standard results imply that for each to > 0, there exists a 
neighborhood / of to on which (11.21) has a unique solution g with q(to) = q(to). To show 
that g = g on J, it suffices to show that g(t) = P{B{t) < q(t)) satisfies g = a on L Since 
g{to) = a, it suffices to show g' = on /. For this, we simply differentiate 

g{t) = / u{x, t) dx, 
J —00 

which gives 

/•<?(*) I I r^it) 

g'{t) = u{q{t),t)q'{t)+ / dtu{x,t) dx = --dM<l{t),t) + - / dlu{x,t) dx = 0. 

j-00 ^ ^ J -00 

Since to > was arbitrary, this shows that q G C°°(0, 00) and satisfies (II. 2p for all t > 0. 
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To show that q(t) — )■ g(0) as t — > 0, we first show that g(0) < hminft^o Qit)- The proof is 
by contradiction. Suppose not. Then there exists e > and t„, J, such that q(tn) < q{0) —e 
for all n. Hence, 

a = P{B{tn) < q{tn)) < P{B{tn) < g(0) - e) 

P(fi(0) < g(0) -e)< P(5(0) < g(0)) = a. 

It follows that P{B{0) < q{0) —e) = a, but this contradicts the uniqueness of the a-quantile 
of the measure f{x) dx. The proof that limsup^^g q{t) < q{0) is similar. □ 

We next present a proof of the quantile central limit theorem. The proof uses a multi- 
dimensional version of the Lindeberg- Feller theorem, which is stated below as Theorem lA.ll 

Theorem A.l. For each fixed n, let {Xm.n}m=i independent, W^-valued random vectors 
with mean zero and covariance matrix am,n- If 

(V Zlm=i (^m,n cT as u ^ oo , and 

(a) for each 9 eM.'^ and each e > 0, J2m=i ^[\^ ' ^m,npl{|e-x™.„|>e}] —'■0 as n ^ oo, 

then Sn = Xi^n + ■ ■ ■ + Xn^n =^ N, where N is multi-normal with mean and covariance a. 

Proof of Theorem 12.11 Given x,y E M.'^, we shall write x < ?/ if x{i) < y{i) for all i. 
Fix X G and for each n, m G N, 1 < m < n, define the random vector Ym^n ^ by 

^m,n(j) = n (l{X,„(i)<n-i/2x(i)+<?(i)} - Pn{j)) , 

where PniJ) = + q{j))- Then for each fixed n G N, Fm,n, 1 < m < n, are 

independent, EYm^n = 0, and 

(i) Em=i E[Ym,n{'^Wni,n{j)] ^ Gij{q{i), q{j)) - a{i)a{j) as n ^ oo, 

(ii) for each ^ G M'^ and e > 0, Zl=i E[\e ■ l"m,„,|'l{|e.y„,„|>.}] -> as n ^ oo. 
Part (i) follows since 

n 

Y,E[Ym,nmmAj)] 
m=l 

n 

= n-' 5^[P(X„(^) < n-'/^x{^ + q{z),XmU) < n-'^'xij) + qU)) - Pni^PnU)] 

m=l 

= P{X(t)<n''/'xit) + q(t),X{j)<n-'/^xij) + q{j))-p^{t)PnU), 

and part (ii) follows since 1^ ■ Ym,n\ < n-i/2^max(|^(l)|, . . . , \9{d)\), and therefore P{\9 ■ 
Ym,n\ > e) = for sufficiently large n. 

Thus, by Theorem I A. 11 Sn. = Yin + ■ ■ ■ + Y^n =^ where TV is multinormal with mean 
and covariance 

E[iV(z)iV(j)] = G,,(g(z),g(j))-a(Oa(j)- (A.l) 
Now, observe that the following are equivalent: 
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(a) n^/^(X«,(„),„ -q) <x, 

(b) < n-i/2x(j) + for all j, 

(C) Em=l l{X™0)<n-l/2x{i)+q{i)} > l^ij^, ]) ioi all J, 

(d) ri-^/2^;;^i(l{^^(^.)<„-i/2^.(j)+^(j)l > n-^l^{K{n,]) - np^U)) for all j. 

Thus, if a„, G M"' is defined by a„(j) = n~^^'^{K{n, j) — npn{j)), then P(?7,^/^(Xk(„).„ — g) < 
= P{Sn > fln)- Note that 

so that a„ — 7- a G M'', where a(j) = — Therefore, 

P(ni/2(X,(„).„ _ g) < a;) ^ P(iV > a) = P(iV < -a) = P{N < x). 



where N is the random vector defined by N{j) = N[j)/^j[q{j)). Comparing (1A.1I) and 



2.ip . this completes the proof. □ 

At certain points in Section El we use the fact that P„ has well-behaved tail probabilities. 
The precise formulation of this fact is given below. 

Proposition A. 2. Fix T > 0. For each r > 0, there exist constants C and uq such that 

P{\Fn{t)\ > A) < C(A"^ + P(|P(t) - q{t)\' > A)), 
for all n > uq, <t <T, and A > 0. In particular, sup„ P(|P„(0)| > A) as A —t- oo. 

Remark A. 3. Since B is Gaussian, Proposition lA.2l in fact shows that P(|P„(t)| > A) < C\~^. 
Indeed, 

P(|P(t) - > A) < P(|P(t)| + M> A^/^) = 2P(P(t) < -A^/' + M), 

where M = supo<t<T If A > 4M^, then 

P(|P(t) - > A) < 2P(P(t) < -A'/V2) = 2$(-t-^/2y/2/2) < 2^{-T-'^^X^/y2), 

where $ is the distribution function of the standard normal, which satisfies $(— x) < CrX~'^ 
for all r > 0. Hence, P(|P(t) - > A) < 2C2,.T-^-2-2^A^' for A sufficiently large. 

Proof of Proposition IA.21 Fix T > and r > 0. Since u{q{0),0) > 0, there exists 
e G (0, 1) such that 

m = inf{M(a:,t) : \x - q{t)\ < 6,0 <t <T} > 0. 

Since j{n)/n = a + o(n~^/^), we may choose rio > 2m~^ such that \j{n)/n — a\ < n~^^'^em/2, 
for all n > Uq. Let n > uq, t ^ [0,T], and A > be arbitrary. Note that we may assume 
without loss of generality that A > 1 + 2m^^/^ and r > 2. 
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Let us first assume n ^/^A < e. We begin with 

P(F„(t) < -A) = p(j2hBAt)«iit)-n--^'M > 

V j=i J 

= P\n~^^'y\{B,{t)<q{t)-n~^l^\} -Pn) > — PnJ, (A.2) 

where Pn = Pn(A) = P{B(t) < q(t) — n^^/'^X). By the mean value theorem, a—pn > n~^/^Am. 
Since \j{n)/n — a\ < n~^/^m/2 < n~^/'^\m/2^ we have 

P(F„(t) < -A) < p(^"'X^(l{B,(i)<,(i)-„-i/2A} - Pn) > n-i/2^m/2Y 

V i=i ^ 

By Chebyshev's inequahty, 

P{Fn{t) < -A) < Cn-"/2^-'^E 

By Burkholder's inequahty (see, for example. Theorem 6.3.10 in |16]). 



X^(l{Bj(i)<g(i)-n-i/2A} - Pn] 



n r/2 
P(F„(t) < -A) < Cn-^/^A-^P J2 |l{B,W<,(t)~n-Vn} -PnP . 

i=i 

Finally, by Jensen's inequality, 

P(F„(t) < -A) < CA-'"P|l|^(,)<,(,)_„-v.A} -Pnr- 

Note that 

= p„(l-p„)((l-p„)^-i+p;-^)<2p„. 

Hence, P(P„,(t) < -A) < CX~'pn{X) < CA'^ 
Similarly, using the mean value theorem, 

P(P„(t) > A) = ^(X^ l{B,W>g(t)+n-i/2A} > ^ - JH + l) 

< PU-^ 2^(l{i?,W>gW+n-i/2A} - Pn) > 1 - ^ - P„ j 

< PU"^^(l{B,(i)>g(i)+n-i/2A} - Pn) > ri"^/2Am/2 j . 

where p„ = p„(A) = P{B(t) > q(t) + n~^/^A). By Chebyshev, Burkholder, and Jensen, 
P{Fn{t) > A) < CA^''p„(A) < CA-^ This completes the proof in the case n-^/'^X < e. 
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We now assume n~^/^A > e. In this case, a — p„ > P{q{t) — e < B(t) < q{t)) > em. 
Hence, j{n)/n —pn > em/2. Thus, starting from (]A.2|) . Chebyshev, Burkholder, and Jensen 
imply 



P(F„(t) < -A) < Cn-'^E 



X^(l{Bj(t)<g(t)-n-l/2A} - Pn) 



2r 



Similarly, P(F„(t) > A) < Cn-^p„(A). Hence, 

P{\Fn{t)\ > A) < Cn"^(p„(A) +p„(A)) = Cn~''P{\B{t) - q{t)\ > n~'/^X). 
If n > A, then P(|F„(t)| > A) < CX~^', and we are done. If n < A, then 

P{\Fn{t)\ > A) < CP{\B{t) - q{t)\ > X'/') = CPmt) - q{t)\' > A), 
and this completes the proof. 



□ 



In the proof of Theorem 15. we need the following formula for the distribution of the 
j-th order statistic. 

Lemma A. 4. If {Xn} is a sequence of iid copies of X, then for all a; G M, 
where = P{X < x) and $ = 1 — $. 

Proof. Recall that $ is continuous. Hence, if g is absolutely continuous on [0, 1], then 



J — oo 

(See Exercise 3.36(b) in [6], for example.) We therefore have 



(n \ n y n 

i=l ' k=j ^ i=l 



E 

k=j 



^ k=j 
\k'^^„,\n—k—l 



n 



,n—k 



{k^{yf-'^{yr~' - (n - k)^{yf^{yr-'-') d^{y). 

J J ~oo 

We can rewrite this as 
P{Xr.n <x) = £^ [gfcQ$(l/)'=-^$(l/)"-'= 



k=j+i 



n 



k - 1 



^y)'~'<l>iyr-' 



d^y). 
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n 
k-l 



k 



Observing that 

{n-k + 1) 

completes the proof. 

The following lemma is needed in the proof of Theorem | 
Lemma A. 5. For each r > 1, there exist constants C > and uq E N such that 



□ 



E 



i=l 



2r 



< C{{npY y (np)), 



(A.3) 



for all n > no and all p G [0, 1]. Note that C depends only on r, and not on n or p. 



Remark A. 6. If r < 1, then Lemma [A. 5 1 together with Jensen's inequality, imply that the 
expectation in (1A.3|) is bounded above by C{npY . 

Proof of Lemma IA.51 First assume that r G N. Let ^ = l{t/j<p} —p and ip{t) = i?[e**^]. 
Then 



E 



i=l 



2r 



t=0 



By Faa di Bruno's formula (see [9], for example), if n > 2r, then 



dt^ 



b,\b2\---b2rHn-k)\ 



1! 



2! 



(2r)! 



b2r 



where the sum is over all different solutions in nonnegative integers 6i, . . . , 62r- of 6i + 262 + 
363 + ■ ■ ■ + 2r62r = 2r, and k = bi + ■ ■ ■ + b2r- Let us take t = and observe that, since 
^'{0) = 0, every nonzero summand has 61 = 0, which implies 



, 262 + 263 + 264 + ■ ■ ■ + 262r / 262 + 363 + 

k = < 



2rh 



2r 



r. 



Also, (^(0) = 1 and |(^(^')(0)| = \Ee.^\ < 2p. Therefore 



dt^' 



t=0 



k=l 



k=l 



Taking into account the two possibilities, np < 1 and up > 1, gives us (]A.3p . 

Now consider the general case, r G [l,oo). First assume np > 1. Choose an integer 
k > r. Then (1A.3I) . together with Jensen's inequality, give 



E 



n 2r 



< E 



n 2k^ 



r/k 



< C{npY. 



i=l 
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Next assume np < 1. Choose positive integers k,i such that i <r < k. Then 



E 



i=l 



2r 



1](1{C/,<P} -P) <E 5^(l{t/,<p} - P) 



i=l 



21 



E 



i=l 



2k 



< Ci{np) + C2{np) = C{np), 



which completes the proof. 

Finally, we present the proof of Lemma 
tightness in the large gap regime. 



□ 



which establishes the needed estimates for 



Proof of Lemma 16.31 Fix T > 0, A G (0, 1/2), and p > 2. Without loss of generality, 
assume s < t and define 5 = t — s. Note that it is sufficient to prove that there exists a 
constant 6o > such that the lemma holds whenever 5 < ^o- 

Let M = supo<t<r and define = (M"^ A 2)/4. Note that, for a.e. u e fl, 

#{j : Bj{s,uj) < Qn{s,u)} = j{n). Hence, if Bj{t,uj) < Bj{s,uj) + rT^I'^e + g(t) — for 
all j, then #{j : Bj{t, oj) < Qnis, cu) + n'^^'^e + q{t) — q{s)} > j{n). In other words, up to a 
set of measure zero, 

n 

f]{B,{t) < B,{s) + n-^l\ + q(t) - C {Q^it) < g„(s) + n-'/^e + q{t) - q{s)}. 

i=i 

We therefore have 

P{E^it) - >e) = PiQnit) - Qnis) > n-^/^ + g(t) - 

< p((j{BAt) - B,{s) > n-'l\ + q{t) - q{s)}\ 

< nP[B{t) - B{s) > n-^l'^e + q{t) - 

< nP{B{t) - B{s) > n-^l'^e - M\t - s\). 

Recalling $, defined below (16. 6p . this gives 

P(F„(t) - >e)< _ Mb'^''). 

Note that Mb^l'^ < MSl''^ < 1/2 < {1/2)5-^ and n~^/^e5^^l'^ > 6~^. Thus, using 
< CrX~^ with r = A"^(l +p/4), we have 

P{Fn{t) - Fn{s) >e)< nCrrS'^ < e^6~^+^Cr2'6^-^P/^ = Cp,Ae^6'^\t - s\p/\ 

which is in fact a sharper bound than necessary. (The second inequality above follows from 
the assumption in the statement of the lemma that n~^^'^e > \t — s\^^'^~^.) The bound on 
P{Fn{t) — Fn{s) < —e) is obtained similarly. □ 
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